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Abstract 

Exclusive p° electroproduction at HERA has been studied with the ZEUS de- 
tector using 120 pb _1 of integrated luminosity collected during 1996-2000. The 
analysis was carried out in the kinematic range of photon virtuality 2 < Q 2 < 160 
GeV 2 , and j*p centre-of-mass energy 32 < W < 180 GeV. The results include 
the Q 2 and W dependence of the YP ~~ * P°P cross section and the distribution 
of the squared-four-momentum transfer to the proton. The helicity analysis of 
the decay-matrix elements of the p° was used to study the ratio of the j*p cross 
section for longitudinal and transverse photon as a function of Q 2 and W. Fi- 
nally, an effective Pomeron trajectory was extracted. The results are compared 
to various theoretical predictions. 
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1 Introduction 



Two of the most surprising aspects of high-energy deep inelastic scattering (DIS) observed 
at the HERA ep collider have been the sharp rise of the proton structure function, F2, 
with decreasing value of Bjorken x and the abundance of events with a large rapidity gap 
in the hadronic final state [1]. The latter are identified as due to diffraction in the deep 
inelastic regime. A contribution to the diffractive cross section arises from the exclusive 
production of vector mesons (VM). 

High-energy exclusive VM production in DIS has been postulated to proceed through 
two-gluon exchange [2,3], once the scale, usually taken as the virtuality Q 2 of the ex- 
changed photon, is large enough for perturbative Quantum Chromodynamics (pQCD) to 
be applicable. The gluons in the proton, which lie at the origin of the sharp increase of 
F 2 , are also expected to cause the VM cross section to increase with increasing photon 
proton centre-of-mass energy, W, with the rate of increase growing with Q 2 . Moreover, 
the effective size of the virtual photon decreases with increasing Q 2 , leading to a flatter 
distribution in t, the four-momentum-transfer squared at the proton vertex. All these 
features, with varying levels of significance, have been observed at HERA [4-10] in the 
exclusive production of p°, u>, 0, and J/tp mesons. 

This paper reports on an extensive study of the properties of exclusive p°-meson produc- 
tion, 

7*p -> p°p, 

based on a high statistics data sample collected with the ZEUS detector during the period 
1996-2000, corresponding to an integrated luminosity of about 120 pb _1 . 

2 Theoretical background 

Calculations of the VM production cross section in DIS require knowledge of the qq wave- 
function of the virtual photon, specified by QED and which depends on the polarisation of 
the virtual photon. For longitudinally polarised photons, 7£, qq pairs of small transverse 
size dominate [3]. The opposite holds for transversely polarised photons, 7^, where qq 
configurations with large transverse size dominate. The favourable feature of exclusive 
VM production is that, at high Q 2 , the longitudinal component of the virtual photon is 
dominant. The interaction cross section in this case can be fully calculated in pQCD [11], 
with two-gluon exchange as the leading process in the high-energy regime. For heavy 
vector mesons, such as the J /if) or the T, perturbative calculations apply even at Q 2 = 0, 
as the smallness of the qq dipole originating from the photon is guaranteed by the mass 
of the quarks. 
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Irrespective of particular calculations [12], in the region dominated by perturbative QCD 
the following features are predicted: 

• the total 7*p — > Vp cross section, <r 7 * p , exhibits a steep rise with W, which can be 
parameterised as a ~ W s , with S increasing with Q 2 ; 

• the Q 2 dependence of the cross-section, which for a longitudinally polarised photon is 
expected to behave as Q~ 6 , is moderated to become Q~ 4 by the rapid increase of the 
gluon density with Q 2 ; 

• the distribution of t becomes universal, with little or no dependence on W or Q 2 ; 

• breaking of the s-channel helicity conservation (SCHC) is expected. 

In the region where perturbative calculations are applicable, exclusive vector-meson pro- 
duction could become a complementary source of information on the gluon content of the 
proton. At present, the following theoretical uncertainties have been identified: 

• the calculation of e>"(7*p — > Vp) involves the generalised parton distributions [13,14], 
which are not well tested; in addition [15], it involves gluon densities outside the range 
constrained by global QCD analyses of parton densities; 

• higher-order corrections have not been fully calculated [16]; therefore the overall nor- 
malisation is uncertain and the scale at which the gluons are probed is not known; 

• the rapid rise of <t 7 * p with W implies a non-zero real part of the scattering amplitude, 
which is not known; 

• the wave-functions of the vector mesons are not fully known. 

In spite of all these problems, precise measurements of differential cross sections sepa- 
rated into longitudinal and transverse components [17], should help to resolve the above 
theoretical uncertainties. 

It is important in these studies to establish a region of phase space where hard interactions 
dominate over the non-perturbative soft component. If the relative transverse momentum 
of the qq pair is small, the colour dipole is large and perturbative calculations do not apply. 
In this case the interaction looks similar to hadron-hadron elastic scattering, described 
by soft Pomeron exchange as in Regge phenomenology [18]. 

The parameters of the soft Pomeron are known from measurements of total cross sections 
for hadron-hadron interactions and elastic proton-proton measurements. It is usually 
assumed that the Pomeron trajectory is linear in t: 

a P (t) = a P (0) + a'jp t . (1) 

The parameter a P (0) determines the energy behaviour of the total cross section, 

a tot ~ (W*)«*(o)-i 
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and a'jp describes the increase of the slope b of the t distribution with increasing W. The 
value of a'p is inversely proportional to the square of the typical transverse momenta 
participating in the exchanged trajectory. A large value of a' p suggests the presence of 
low transverse momenta typical of soft interactions. The accepted values of a P (0) [19] 
and a'p [20] are 



« P (0) = 1.096 ±0.003 
a' P = 0.25 GeV" 2 . 

The non-universality of a P (0) has been established in inclusive DIS, where the slope of 
the 7*p total cross section with W has a pronounced Q 2 dependence [21]. The value of 
a' F can be determined from exclusive VM production at HERA via the W dependence of 
the exponential b slope of the t distribution for fixed values of W, where b is expected to 
behave as 

b(W)=b + 4a P \n^-, 
Wo 

where b and Wq are free parameters. The value of a' p can also be derived from the W 
dependence of da/dt at fixed t, 

^-(W) = F{t)W^ 2a ^- 2 \ (2) 

LLC 

where F(t) is an arbitrary function. This approach has the advantage that no assumption 
needs to be made about the t dependence. The first indications from measurements of 
a F (t) in exclusive J/ip photoproduction [8,22] are that a P (0) is larger and a' p is smaller 
than those of the above soft Pomeron trajectory. 



3 Experimental set-up 

The present measurement is based on data taken with the ZEUS detector during two 
running periods of the HERA ep collider. During 1996-1997, protons with energy 820 
GeV collided with 27.5 GeV positrons, while during 1998-2000, 920 GeV protons collided 
with 27.5 GeV electrons or positrons. The sample used for this study corresponds to an 
integrated luminosity of 118.9 pb" 1 , consisting of 37.2 pb" 1 e + p sample from 1996-1997 
and 81.7 pb^ 1 from the 1998-2000 sample (16.7 pb^ 1 e~ and 65.0 pb" 1 e+jj 

A detailed description of the ZEUS detector can be found elsewhere [23,24]. A brief 
outline of the components that are most relevant for this analysis is given below. 

1 From now on, the word "electron" will be used as a generic term for both electrons and positrons. 
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Charged particles are tracked in the central tracking detector (CTD) [25-27]. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar-angl^l region 15° < 6 < 164°. The CTD operates in a magnetic field of 1.43 T 
provided by a thin solenoid. The transverse-momentum resolution for full-length tracks 
is a{p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [28-31] covers 99.7% of the 
total solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) 
and the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section (EMC) and either one (in RCAL) or two 
(in BCAL and FCAL) hadronic sections. The CAL energy resolutions, as measured under 
test-beam conditions, are a(E)/E = 0.18 /\/~E for electrons and a(E)/E = 0.35/v^ for 
hadrons, with E in GeV. 

The position of the scattered electron was determined by combining information from the 
CAL, the small-angle rear tracking detector [32] and the hadron-electron separator [33]. 

In 1998, the forward plug calorimeter (FPC) [34] was installed in the 20x20 cm 2 beam 
hole of the FCAL with a small hole of radius 3.15 cm in the centre to accommodate the 
beam pipe. The FPC increased the forward calorimeter coverage by about one unit in 
pseudorapidity to rj<5. 

The leading-proton spectrometer (LPS) [35] detected positively charged particles scat- 
tered at small angles and carrying a substantial fraction, xl, of the incoming proton 
momentum; these particles remained in the beam-pipe and their trajectories were mea- 
sured by a system of silicon microstrip detectors, located between 23.8 m and 90.0 m 
from the interaction point. The particle deflections induced by the magnets of the proton 
beam-line allowed a momentum analysis of the scattered proton. 

During the 1996-1997 data taking, a proton- remnant tagger (PRT1) was used to tag 
events in which the proton dissociates. It consisted of two layers of scintillation counters 
perpendicular to the beam at Z = 5.15 m. The two layers were separated by a 2 mm-thick 
lead absorber. The pseudorapidity range covered by the PRT1 was 4.3 < rj < 5.8. 

The luminosity was measured from the rate of the bremsstrahlung process ep —>■ ejp. 
The photon was measured in a lead-scintillator calorimeter [36-38] placed in the HERA 
tunnel at Z = —107 m. 



2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton direction, referred to as the "forward direction" , and the X axis pointing left towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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4 Data selection and reconstruction 



The following kinematic variables are used to describe exclusive p production and its 
subsequent decay into a 7r + 7r~ pair: 

• the four-momenta of the incident electron (&), scattered electron (k'), incident proton 
(P), scattered proton (P') and virtual photon (q); 

• Q 2 = —q 2 = — [k — k') 2 , the negative squared four-momentum of the virtual photon; 

• W 2 = (q + P) 2 , the squared centre-of-mass energy of the photon-proton system; 

• y = (P • q)/(P • k), the fraction of the electron energy transferred to the proton in its 
rest frame; 

• M^,,-, the invariant mass of the two decay pions; 

• t = (P — P') 2 , the squared four-momentum transfer at the proton vertex; 

• three helicity angles, 6h and 4>h ( see Section [9]). 

The kinematic variables were reconstructed using the so-called "constrained" method [10, 
39], which uses the momenta of the decay particles measured in the CTD and the recon- 
structed polar and azimuthal angles of the scattered electron. 

The online event selection required an electron candidate in the CAL, along with the 
detection of at least one and not more than six tracks in the CTD. 

In the offline selection, the following further requirements were imposed: 

• the presence of a scattered electron, with energy in the CAL greater than 10 GeV 
and with an impact point on the face of the RCAL outside a rectangular area of 
26.4 x 16 cm 2 ; 

• E — Pz > 45 GeV, where E — Pz = ^2%i^i — pzj and the summation is over the 
energies and longitudinal momenta of the final-state electron and pions, was imposed. 
This cut excludes events with high energy photons radiated in the initial state; 

• the Z coordinate of the interaction vertex within ±50 cm of the nominal interaction 
point; 

• in addition to the scattered electron, exactly two oppositely charged tracks, each 
associated with the reconstructed vertex, and each having pseudorapidity \r)\ less than 
1.75 and transverse momentum greater than 150 MeV; this excluded regions of low 
reconstruction efficiency and poor momentum resolution in the CTD. These tracks 
were treated in the following analysis as a 7T + n~ pair; 

• events with any energy deposit larger than 300 MeV in the CAL and not associated 
with the pion tracks (so-called 'unmatched islands') were rejected [40-42]. 
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In addition, the following requirements were applied to select kinematic regions of high 
acceptance: 

• the analysis was restricted to the kinematic regions 2 < Q 2 < 80 GeV 2 and 32 < W < 
160 GeV in the 1996-1997 data and 2 < Q 2 < 160 GeV 2 and 32 < W < 180 GeV in 
the 1998-2000 sample; 

• only events in the n + ir~ mass interval 0.65 < < 1.1 GeV and with \t\ < 1 GeV 2 
were taken. The mass interval is slightly narrower than that used previously [10], in 
order to reduce the effect of the background from non-resonant tt + it~ production. In 
the selected M n7T range, the resonant contribution is ~ 100% (see Section [8]). 

The above selection yielded 22,400 events in the 1996-1997 sample and 49,300 events in 
the 1998-2000 sample, giving a total of 71,700 events for this analysis. 



5 Monte Carlo simulation 

The relevant Monte Carlo (MC) generators have been described in detail previously [10]. 
Here their main features are summarised. 

The program Zeusvm [43] interfaced to Heracles4.4 [44] was used. The effective Q 2 , 
W and t dependences of the cross section were parameterised to reproduce the data [42]. 

The decay angular distributions were generated uniformly and the MC events were then 
iteratively reweighted using the results of the present analysis for the 15 combinations of 
matrix elements r^, rf k (see Section [9]). 

The contribution of the proton-dissociative process was studied with the Epsoft [45] 
generator for the 1996-1997 data and with Pythia [46] for the 1998-2000 data. The Q 2 , 
W and t dependences were parameterised to reproduce the control samples in the data. 
The decay angular distributions were generated as in the Zeusvm sample. 

The generated events were processed through the same chain of selection and reconstruc- 
tion procedures as the data, thus accounting for trigger as well as detector acceptance and 
smearing effects. For both MC sets, the number of simulated events after reconstruction 
was about a factor of seven greater than the number of reconstructed data events. 

All measured distributions are well described by the MC simulations. Some examples are 
shown in Fig. [U for the W, Q 2 , t variables, and the three helicity angles, Oh, (frh, and 
and in Fig. [2] for the transverse momentum pt of the pions, for different Q 2 bins. 
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6 Systematics 



The systematic uncertainties of the cross section were evaluated by varying the selection 
cuts and the MC simulation parameters. The following selection cuts were varied: 

• the E — Pz cut was changed within the appropriate resolution of ±3 GeV; 

• the pt of the pion tracks (default 0.15 GeV) was increased to 0.2 GeV; 

• the distance of closest approach of the extrapolated track to the matched island in the 
CAL was changed from 30 cm to 20 cm; 

• the 7r + 7r _ -mass window was changed to 0.65-1.2 GeV; 

• the Z vertex cut was varied by ±10 cm; 

• the rectangular area of the electron impact point on the CAL was increased by 0.5 cm 
in X and Y; 

• the energy of an unmatched island was lowered to 0.25 GeV and then raised to 0.35 
GeV. 

The dependence of the results on the precision with which the MC reproduces the perfor- 
mance of the detector and the data was checked by varying the following inputs within 
their estimated uncertainty: 

• the reconstructed position of the electron was shifted with respect to the MC by 
±1 mm; 

• the electron-position resolution was varied by ±10% in the MC; 

• the W^-dependence in the MC was changed by varying 5 by ±0.03; 

• the exponential t-distribution in the MC was reweighted by changing the nominal 
slope parameter b by ±0.5 GeV" 2 ; 

• the angular distributions in the MC were reweighted assuming SCHC; 

• the Q 2 -distribution in the MC was reweighted by (Q 2 ± M 2 ) k , where k = ±0.05. 

The largest uncertainty of about ±4% originated from the variation of the energy of 
the unmatched islands. All the other checks resulted on average in a 0.5% change in 
the measured cross sections. All the systematic uncertainties were added in quadrature. 
In addition, the cross-section measurements have an overall normalisation uncertainty 
of ±2% due to the luminosity measurement. 
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7 Proton dissociation 



The production of p° mesons may be accompanied by the proton-dissociation process, 
7*p — > p°N. For low masses Mn of the dissociative system N, the hadronisation products 
may remain inside the beam-pipe, leaving no signals in the main detector. The contribu- 
tion of these events to the exclusive p° cross section was estimated from MC generators 
for proton- dissociative processes. 

A class of proton dissociative events for which the final-state particles leave observed 
signals in the surrounding detectors was used to tune the and the t distribution 
in the MC. In the 1998-2000 running period, these events were selected by requiring 
a signal in the FPC detector with energy above 1 GeV. The comparison of the data 
with Pythia expectations for the energy distribution in the FPC is shown in Fig. [3](a). 
The same procedure was repeated with a sample of p° events for which the FPC energy 
was less than 1 GeV and a leading proton was measured in the LPS detector, with 
the fraction of the incoming proton momentum Xl < 0.95. The comparison between 
the xl distribution measured in the data and that expected from Pythia is shown in 
Fig. E](b), where the elastic peak in the data (x^ > 0.95) is also observed. Also shown 
in Fig. Elc-e) is the fraction of proton-dissociative events expected in the selected p° 
sample as a function of Q 2 , W and t. The fraction is at the level of 19%, independent 
of Q 2 and W, but increasing with increasing \t\. The combined use of the FPC and LPS 
methods leads to an estimate of the proton dissociative contribution for \t\ < 1 GeV 2 of 
0.19 ± 0.02(stat.) ± 0.03(syst.). The systematic uncertainty was estimated by varying 
the parameters of the Mm distribution and by changing the FPC cut. 

In the 1996-1997 data-taking period, a similar procedure was applied, after tuning the 
Epsoft MC to reproduce events with hits in the PRT1 or energy deposits in the FCAL. 
The proton-dissociative contribution for \t\ < 1 GeV 2 was determined to be 0.07 ± 0.02 
after rejecting events with hits in the PRT1 or energy deposits in the FCAL. This number 
is consistent with that determined from the LPS and FPC because of the different angular 
coverage of the PRT1. 

After subtraction of the proton-dissociative contribution, a good agreement between the 
cross sections derived from the two data-taking periods was found. For all the quoted cross 
sections integrated over t, the overall normalisation uncertainty due to the subtraction of 
the proton-dissociative contributions was estimated to be ±4% and was not included in the 
systematic uncertainty. The proton-dissociative contribution was statistically subtracted 
in each analysed bin, unless stated otherwise. 
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8 Mass distributions 



The 7r + 7r~-invariant-mass distribution is presented in Fig. |H A clear enhancement in the 
p° region is observed. Background coming from the decay <fi — > K + K~, where the kaons 
are misidentified as pions, is expected [42] in the region M njT < 0.55 GeV. That coming 
from uj events in the decay channel uj — > 7r + 7r~7r°, where the ir° remains undetected, 
contributes [42] in the region M nn < 0.65 GeV. Therefore defining the selected p° events 
to be in the window 0.65 < M nn < 1.1 GeV ensures no background from these two 
channels. 

In order to estimate the non-resonant 7r + 7r~ background under the p°, the Soding param- 
eterisation [47] was fitted to the data, with results shown in the figure. The resulting 
mass and width values are in agreement with those given in the Particle Data Group [48] 
compilation. The integrated non- resonant background is of the order of 1% and is thus 
neglected. 

The n + 7r~ mass distributions in different regions of Q 2 and t are shown in Fig. [5] and 
Fig. O respectively. The shape of the mass distribution changes neither with Q 2 nor with 
t. The results of the fit to the Soding parameterisation are also shown. Note that the 
interference term decreases with Q 2 as expected but is independent of t, indicating that 
the non-exclusive background is negligible. 

9 Angular distributions and decay-matrix density 

The exclusive electroproduction and decay of p° mesons is described, at fixed W, Q 2 , M n7T 
and t, by three helicity angles: $^ is the angle between the p° production plane and the 
electron scattering plane in the 7*p centre-of-mass frame; Oh and <ph are the polar and 
azimuthal angles of the positively charged decay pion in the s-channel helicity frame. In 
this frame, the spin-quantisation axis is defined as the direction opposite to the momen- 
tum of the final-state proton in the p° rest frame. In the 7*p centre-of-mass system, <ph is 
the angle between the decay plane and the p° production plane. The angular distribution 
as a function of these three angles, W (cos Oh, 4>hi ®h), is parameterised by the p° spin- 
density matrix elements, pf k , where i,k = —1,0,1 and by convention a—0, 1,2,4,5,6 for 
an unpolarised charged-lepton beam [49]. The superscript denotes the decomposition of 
the spin-density matrix into contributions from the following photon-polarisation states: 
unpolarised transverse photons (0); linearly polarised transverse photons (1,2); longitudi- 
nally polarised photons (4); and from the interference of the longitudinal and transverse 
amplitudes (5,6). 

The decay angular distribution can be expressed in terms of combinations, r®£ and r" k , of 
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the density matrix elements 




1 + eR 

Ptk 



a 



1,2 



1 + eR' 



r 



< 



ik 



. 1 + eR' 



a 



5,6, 



where e is the ratio of the longitudinal- to transverse-photon fluxes and R = ctl/vt, with 
cl and ot the cross sections for exclusive p° production from longitudinal and transverse 
virtual photons, respectively. In the kinematic range of this analysis, the value of e 
varies between 0.96 and 1 with an average value of 0.996; hence p\ k and p\ k cannot be 
distinguished. 

The Hermitian nature of the spin-density matrix and the requirement of parity conser- 
vation reduces the number of independent parameters to 15 [49]. A 15-parameter fit was 
performed to the data and the obtained results are listed in Table [Hand shown in Fig. [7] as 
a function of Q 2 . The published ZEUS results [50] at lower Q 2 values and the expectations 
of SCHC, when relevant, are also included. The observed Q 2 dependence, expected in 
some calculations [51] and previously reported by HI [52], is driven by the R dependence 
on Q 2 under the assumption of helicity conservation and natural parity exchange. The 
significant deviation of Tq from zero shows that SCHC does not hold [51] as was observed 
previously [50,52]. 

The angular distribution for the decay of the p° meson, integrated over <ph and $>h, reduces 
to 



The element 7qq may be extracted from a one-dimensional fit to the cos 6h distribution. 
The cos 9h distributions, for different Q 2 intervals, are shown in Fig. [HJ together with the 
results of a one-dimensional fit of the form ([3]). The data are well described by the fitted 
parameter at each value of Q 2 . 

10 Cross section 

The measured 7*p cross sections are averaged over intervals listed in the appropriate tables 
and are quoted at fixed values of Q 2 and W. The cross sections are corrected for the mass 
range 0.28 < M wn < 1.5 GeV and integrated over the full t-range, where applicable. 



W(cos0 h ) oc [(1 



ro°o) + (3r o 4 -l)cos 2 ^]. 



(3) 



10 



10.1 t dependence of ct{Yp ~^ P°p) 

The determination of er(7*p — » p°p) as a function of t for W = 90 GeV was performed 
by averaging over 40 < W < 140 GeV. The differential cross-section da/dt{^f*p — » 
is shown in Fig. [9] and listed in Table [2], for different ranges of Q 2 . An exponential form 
proportional to e~ fe '*' was fitted to the data in each range of Q 2 ; the results are shown in 
Fig. [TUJ The exponent b, listed in Table [3], decreases as a function of Q 2 . After including 
the previous results at lower Q 2 [10,53], a sharp decrease of b is observed at low Q 2 ; the 
value of b then levels off at about 5 GeV" 2 . 

A compilation of the value of the slope b for exclusive VM electroproduction, as a function 
of Q 2 + M 2 , is shown in Fig. (TTJ Here M is the mass of the corresponding final state. 
It also includes the exclusive production of a real photon, the deeply virtual Compton 
scattering (DVCS) measurement [54]. When b is plotted as a function of Q 2 + M 2 , the 
trend of b decreasing with increasing scale to an asymptotic value of 5 GeV -2 , seems to 
be a universal property of exclusive processes, as expected in perturbative QCD [2]. 

10.2 Q 2 dependence of cr{^*p —> p°p) 

The determination of <j(^*p — > p°p) as a function of Q 2 for W = 90 GeV was performed by 
averaging over 40 < W < 140 GeV. The results are shown in Fig. [12] with corresponding 
values given in Table HI As expected, a steep decrease of the cross section with Q 2 
is observed. The photoproduction and the low-Q 2 (< 1 GeV 2 ) measurements are also 
shown in the figure. An attempt to fit the Q 2 dependence with a simple propagator term 

with the normalisation and n as free parameters, failed to produce results with an accept- 
able x 2 - The data appear to favour an n value which increases with Q 2 . 

10.3 W dependence of cr{^*p — > p°p) 

The values of the cross section <r(7*p — > p°p) as a function of W, for fixed values of Q 2 , 
are plotted in Fig. [13] and given in Table [5] The cross sections increase with increasing 
W, with the rate of increase growing with increasing Q 2 . 

In order to quantify the rate of growth and its significance, the W dependence for each 
Q 2 value was fitted to the functional form 

a rsj W 5 . 
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The resulting 5 values are presented as a function of Q 2 in Fig. [TH and listed in Table El 
For completeness, the 8 values from lower Q 2 are also included. A clear increase of S with 
Q 2 is observed. Such an increase is expected in pQCD, and reflects the change of the 
\ow-x gluon distribution of the proton with Q 2 . 

To facilitate the comparison, the ZEUS cross-section data as a function of W have been 
replotted in the Q 2 bins used by HI [9]. The results are shown in Fig. [TSJ The agreement 
between the two measurements is reasonable. However, in some Q 2 bins the shape of the 
W dependence is somewhat different. 

A compilation of the value of the slope 5 for exclusive VM electroproduction, as a function 
of Q 2 + M 2 , is shown in Fig. [161 It also includes the DVCS result [54]. When plotted as 
a function of Q 2 + M 2 , the value of 8 and its increase with the scale are similar for all the 
exclusive processes, as expected in perturbative QCD [2]. 



11 R = cfl/cft and r{Jo 

The SCHC hypothesis implies that r\_ x = — Im{r^_ 1 } and Re{r^ } = — Im{r^ }. In this 
case, the ratio R = o^j can be related to the matrix element, 

R = -T^i (4) 
e 1 - r$ 

and thus can be extracted from the $h distribution alone. 

If the SCHC requirement is relaxed, then the relation between R and Tqq is modified, 

1 r 04 _ A 2 
'00 ^ 



R el-(rg*-A= 
with 



O r 04 
z ' 00 



In the kinematic range of the measurements presented in this paper, the non-zero value 
of A implies a correction of ~3% on R up to the highest Q 2 value, where it is ~10%, and 
is neglected. 

Under the assumption that Eq. (j4j) is valid and for values of e studied in this paper, 
< e >=0.996, the matrix element may be interpreted as 

04 / 

where <7 to t = &L + o"t- When the value of is close to one, as is the case for this analysis, 
the error on R becomes large and highly asymmetrical. It is then advantageous to study 
the properties of Tqq itself which carries the same information, rather than R. 
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The Q 2 dependence of for W = 90 GeV, averaged over the range 40 < W < 140 GeV, 
is shown in Fig. [T7| and listed in Table [7] together with the corresponding R values. The 
figure includes three data points at lower Q 2 from previous studies [10,53]. An initial 
steep rise of with Q 2 is observed and above Q 2 ~ 10 GeV 2 , the rise with Q 2 becomes 
milder. At Q 2 = 40 GeV 2 , constitutes about 90% of the total 7*p cross section. 

The comparison of the HI and ZEUS results is presented in Fig. [18] in terms of the ratio 
R. The HI measurements are at W = 75 GeV and those of ZEUS at W = 90 GeV. Given 
the fact that R seems to be independent of W (see below) , both data sets can be directly 
compared. The two measurements are in good agreement. 

The dependence of R on M n7T is presented in Fig. [19] for two Q 2 intervals. The value of 
R falls rapidly with M nn above the central p° mass value. Although a change of R with 
M n7r was anticipated to be ~ 10% [55], the effect seen in the data is much stronger. The 
effect remains strong also at higher Q 2 , contrary to expectations [55]. Once averaged over 
the p° mass region, the main contribution to R comes from the central p° mass value. 

The W dependence of Tqq, for different values of Q 2 , is shown in Fig. [20] and listed in 
Table [8] Within the measurement uncertainties, r^j is independent of W, for all Q 2 
values. This implies that the W behaviour of <jl is the same as that of cry, a result 
which is somewhat surprising. The qq configurations in the wave function of 7^ have 
typically a small transverse size, while the configurations contributing to 7^ may have 
large transverse size. The contribution to cry of large-size qq configurations, which are 
more hadron-like, is expected to lead to a shallower W dependence than in case of a^. 
Thus, the result presented in Fig. [20] suggests that the large-size configurations of the 
transversely polarised photon are suppressed. 

The above conclusion can also explain the behaviour of 7*qq clS cL function of t, shown in 
Fig. EH and presented in Table M for two Q 2 values. Different sizes of interacting objects 
imply different t distributions, in particular a steeper dar/dt compared to dai/dt. This 
turns out not to be the case. In both Q 2 ranges, Tqq is independent of t, reinforcing the 
earlier conclusion about the suppression of the large-size configurations in the transversely 
polarised photon. 

12 Effective Pomeron trajectory 

An effective Pomeron trajectory can be determined from exclusive p° electroproduction 
by using Eq. Since the W dependence of the proton-dissociative contribution was es- 
tablished to be the same as the exclusive p° sample, no subtraction for proton-dissociative 
events was performed. 
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A study of the W dependence of the differential da/dt cross section at fixed t results 
in values of a P (t), listed in Table [TU1 and displayed in Fig. [22], for Q 2 = 3 GeV 2 (upper 
plot) and 10 GeV 2 (lower plot). A linear fit of the form of Eq. ([TTh shown in the figures, 
yields values of a P (0) and a' p shown in Fig. [23] and listed in Table [TTj The value 
of ctp(0) increases slightly with Q 2 , while the value of a' p is Q 2 independent, within the 
measurement uncertainties. Its value tends to be lower than that of the soft Pomeron [56]. 

An alternative way of measuring the slope of the Pomeron trajectory is to study the 
W dependence of the b slope, for fixed Q 2 values. Figure EH displays the values of b 
as a function of W for two Q 2 intervals (see also Table [T2l) . The curves are a result 
of fitting the data to the expression b = b + 4a^ln(W/Wo)- The resulting slopes of 
the trajectory are a' P = 0.15 ± 0.04 (stat.) ±g;gg (syst.) for < Q 2 > = 3.5 GeV 2 and 
a' p = 0.04 ± 0.06 (stat.) (syst.) for < Q 2 > = 11 GeV 2 . These results are consistent 
with those presented in Table [TTJ. 

13 Comparison to models 

In this section, predictions from several pQCD-inspired models are compared to the mea- 
surements. 

13.1 The models 

All models are based on the dipole representation of the virtual photon, in which the 
photon first fluctuates into a qq pair (the colour dipole), which then interacts with the 
proton to produce the p°. The ingredients necessary in such calculations are the virtual- 
photon wave-function, the dipole-proton cross section, and the p° wave-function. The 
photon wave-function is known from QED. The models differ in the treatment of the 
dipole-proton cross section and the assumed p° wave-function. 

The models of Frankfurt, Koepf and Strikman (FKS) [57,58] and of Martin, Ryskin and 
Teubner (MRT) [59, 60] are based on two-gluon exchange as the dominant mechanism 
for the dipole-proton interaction. The gluon distributions are derived from inclusive 
measurements of the proton structure function. In the FKS model, a three-dimensional 
Gaussian is assumed for the p° wave-function, while MRT use parton-hadron duality and 
normalise the calculations to the data. For the comparison with the present measurements 
the MRST99 [61] and CTEQ6.5M [62] parameterisations for the gluon density were used. 

Kowalski, Motyka and Watt (KMW) [63] use an improved version of the saturation 
model [64, 65] , with an explicit dependence on the impact parameter and DGLAP [66-69] 
evolution in Q 2 , introduced through the unintegrated gluon distribution [70]. Forshaw, 
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Sandapen and Shaw (FSS) [71] model the dipole-proton interaction through the exchange 
of a soft [56] and a hard [72] Pomeron, with (Sat) and without (Nosat) saturation, and 
use the DGKP and Gaussian p° wave-functions. In the model of Dosch and Ferreira 
(DF) [73], the dipole cross section is calculated using Wilson loops, making use of the 
stochastic vacuum model for the non-perturbative QCD contribution. 

While the calculations based on two-gluon exchange are limited to relatively high-Q 2 
values (typically ~ 4 GeV 2 ), those based on modelling the dipole cross section incorporate 
both the perturbative and non-perturbative aspects of p° production. 



13.2 Comparison with data 

The different predictions discussed above are compared to the Q 2 dependence of the cross 
section in Fig. None of the models gives a good description of the data over the 
full kinematic range of the measurement. The FSS model with the three-dimensional 
Gaussian p° wave-function describes the low-Q 2 data very well, while the KMW and DF 
models describe the Q 2 > 1 GeV 2 region well. 

The various predictions are also compared with the W dependence of the cross section, for 
different Q 2 values, in Fig. [2HJ Here again, none of the models reproduces the magnitude 
of the cross section measurements. The closest to the data, in shape and magnitude, 
are the MRT model with the CTEQ6.5M parametrisation of the gluon distribution in 
the proton and the KMW model. The KMW model gives a good description of the Q 2 
dependence of 5, as shown in Fig. [271 

The dependence of b on Q 2 is given only in the FKS and the KMW models as shown in 
Fig. [28j The FKS expectations are somewhat closer to the data. 

The expected Q 2 dependence of Tqq is compared to the measurements in Fig. [291 The MRT 
prediction, using the CTEQ6.5M gluon density, is the only prediction which describes the 
data in the whole Q 2 range. While all the models exhibit a mild dependence of on W, 
consistent with the data as shown in Figs. [3D] and EH none of them reproduces correctly 
the magnitude of in all the Q 2 bins. 

In summary, none of the models considered above is able to describe all the features of 
the data presented in this paper. The high precision of the measurements can be used to 
refine models for exclusive p° electroproduction. 
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14 Summary and Conclusions 



Exclusive p° electroproduction has been studied by ZEUS at HERA in the range 2 < Q 2 < 
160 GeV 2 and 32 < W < 180 GeV with a high statistics sample. The Q 2 dependence of 
the 7*p — > p°p cross section is a steeply falling function of Q 2 . The cross section rises with 
W and its logarithmic derivative in W increases with increasing Q 2 . The exponential slope 
of the t distribution decreases with increasing Q 2 and levels off at about b = 5 GeV -2 . 
The decay angular distributions of the p° indicate s-channel helicity breaking. The ratio 
of cross sections induced by longitudinally and transversely polarised virtual photons 
increases with Q 2 , but is independent of W and of \t\, suggesting suppression of large-size 
configurations of the transversely polarised photon. The effective Pomeron trajectory, 
averaged over the full Q 2 range, has a larger intercept and a smaller slope than those 
extracted from soft interactions. All these features are compatible with expectations of 
perturbative QCD. However, none of the available models which have been compared to 
the measurements is able to reproduce all the features of the data. 
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Element 


2 < g 2 < 3 GeV 2 


3 < Q 2 < 4 GeV 2 


4 < g 2 < 6 GeV 2 


6 < g 2 < 10 GeV 2 


10 < g 2 < 100 GeV 2 


r 04 
'00 


0.590±0.006lg;gjg 


0.659±0.008in?5 


0.725±0.008ig;8£f 


0.752±0.008lg;85J 


0.814±0.010tjj;{gj] 


Re(r? 4 ) 


0.024±0.005±g;gg| 


0.025±0.007lg;oo9 


0.007±0.007±8;g?| 


0.014±0.007i^ 


0.014±0.009±g;g5? 


r 04 
' 1-1 


-o.oo9±o.oortg:8w 


-0.010±0.008±g;g?6 


0.000±0.007+^^ 


-0.016±0.007lg;gJt 


-0.001±0.010j£S£ 




-0.008d=0.007±g;g?| 


-0.023±0.008lH?6 


-0.015±0.008+°;°}° 


-0.032±0.008+°;°^ 


-0.002±0.01 l±agio 


r oo 


-0.037±0.019lffi 


-0.014±0.026+°;°^ 


0.020±0.028+°;j£ 2 


0.019±0.030+°;°^ 


-0.018±0.042±°;°g 


Re(r 1 1 ) 


-0.032±0.007lg;gJ| 


-0.023±0.010±g;gg| 


-0.016±0.009±g;gi| 


-0.006±0.011+°;°^ 


-0.042±0.016+°;° 2 ^ 


r l-i 


0.195±0.009j:8;8i9 


0.151±0.01llg;g^ 


o.i2i±o.oiilg:on 


0.095±0.01llg:rag 


0.100±0.016±g;g§l 


Im(r 2 ) 


0.040±0.007±H 


0.024±0.009lg;85o 


0.029±0.009ig;8ii 


0.031±0.009i^2 


0.026±0.01518:mI 


In^r 2 ^) 


-0.186±0.009+°;°^ 


-0.148d=0.011±g;gi§ 


-0.124±0.012+°;°^ 


-0.107±0.011+°;°^ 


-0.052±0.016lg;Si2 




o.oi8±o.oo3±g;SS5 


0.018±0.004lg;g^ 


0.007±0.003l8;g$ 


0.018±0.004+°;°^ 


0.004±0.005+°;°^ 


'oo 


0.085±0.009±g:S?5 


0.089±0.013lg:Sw 


0.106±0.013+°;°}° 


0.093±0.013+°;°^ 


0.168±0.018+°;°^ 


Re(rf ) 


0.167±0.003lg;g^ 


0.164±0.004lg:ooi 


0.143±0.005±g;g^ 


0.132±0.005lg;gg| 


0.1 10±0.007±g;SSs 


r\-i 


0.000±0.005±g;gg| 


-0.006±0.006lg:8o6 


0.001±0.005±g:gg| 


0.000±0.006lg;gJ| 


O.OOliO.OOTlgiSa 


lm(r%) 


-0.157±0.003tHo4 


-0.147±0.004lHo7 


-0.145±0.004+°;°^ 


-0.135±0.004+°;°^ 


-0.125±0.006j£ora 


Im(r^_ 1 ) 


0.010d=0.005±g;g?| 


-0.005=b0.005±g;gg| 


-o.ooi±o.oo5ig:S?7 


0.008±0.005+°;°^ 


-0.002±0.007+°;^ 



Table 1: Spin density matrix elements for electroproduction of p , for different 
intervals of Q 2 . The first uncertainty is statistical, the second systematic. 



O 2 bin 

(GeV 2 ) 


O 2 
(C^eV 2 ) 

\ VJT v V J 


\ L \ 

(CrPV 2 ) 

1 VJC V 1 


da/dt 
(nb/GeV 2 ) stat. 


D V D u . 


2-4 


2.7 


0.05 




-I-4Q ^ 


+117.3 
-155.3 


9—4 


9 7 


n 1 ^ 


1 98/I 9 


-1-39 8 


+65.4 
-87.7 


9—4 


9 7 


n 9Q 


4^n 7 


-1-1 3 K 
ZEIO.O 


+30.8 
-39.1 


9—4 
9—4 


9 7 
9 7 


53 
83 


1 97 R 
IZ ( .o 

98 1 

zo. i 


ZLU . Z 

-I-3 3 
ZtO.O 


+17.2 
-17.0 
+10.3 
-5.1 


4—6 5 


5 


05 


849 7 




+33.3 
-40.5 


4— ^ 


^ 
^ 


1 ^ 

U.l o 

9Q 


41^8 
loy.o 


-4-1 ^ 4 
4-7 n 

31 1 .U 


+18.9 
-26.1 
+10.6 
-13.8 


4— 5 


5 


^ 


43 7 


-1-3 9 


+5.7 
-5.8 


4—6 5 


5 


83 


1 9 t; 


-1-1 8 

Ztl.O 


+2.2 
-2.2 


6 5— 1 n 


7 8 


05 


4 


-1-1 n 8 

ZLlU.O 


+15.4 
-15.0 


^— 1 D 


7 8 


1 ^ 


156.2 


±7.4 


+5.3 
-13.3 


U. U 1U 


7 8 


9Q 


fi7 3 


+3 3 
ZlO.O 


+4.9 
-4.7 


6.5-10 


7.8 


0.53 


99 1 


ZLl.O 


+2.3 
-3.1 


6 5—1 


7 8 


83 


^ n3 
o.uo 


-4-fl Q4 
zizu.y^ 


+1.48 
-0.92 


1 0—1 ^ 

-LU -LU 


1 1 Q 


0^ 


1 1 s n 

1 lO.U 


_l_ c; pi 


+5.5 
-5.7 


1 0—1 K 


1 1 Q 


1 F5 


7n 9 


-I-3 Q 
zto.y 


+5.2 
-3.6 


1 0—1 5 


1 1 Q 


9Q 


26.8 


±1.7 


+1.7 
-2.6 


1 0—1 5 


1 1 Q 


53 


8.40 


±0.76 


+0.97 
-1.36 


10-15 
15-30 


11.9 
19.7 


0.83 
0.05 


9 fi7 
3Q R 

oy.o 


—l—O c.1 
-1-9 9 


+0.48 
-0.52 
+1.7 
-3.3 


1 ^— 3n 


1 Q 7 


n i ^ 


20 4 


±1 5 


+1.9 
-1.4 


1 ^— 3n 


1 Q 7 


n 9Q 


9.12 


±0.71 


+0.59 
-0.94 


15-30 


19.7 


0.53 


2.73 


±0.31 


+0.39 
-0.38 


15-30 


19.7 


0.83 


0.84 


±0.19 


+0.19 
-0.30 


30-80 
30-80 
30-80 


41.0 
41.0 
41.0 


0.05 
0.15 
0.29 


5.44 
2.28 
1.40 


±0.83 
±0.50 
±0.26 


+0.76 
-0.80 
+0.37 
-0.54 
+0.26 
-0.35 


30-80 


41.0 


0.53 


0.42 


±0.11 


+0.07 
-0.11 


30-80 


41.0 


0.83 


0.15 


±0.07 


+0.06 
-0.07 



Table 2: The differential cross-section da/dt for the reaction 7*p — > p°p for 
different Q 2 intervals. The first column gives the Q 2 bin, while the second column 
gives the Q 2 value at which the cross section is quoted. The normalisation uncer- 
tainty due to luminosity (±2%) and proton-dissociative background (±4%); ^ s n °t 
included. 
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Q 2 bin (GeV 2 ) 


Q 2 (GeV 2 ) 


b (GeV" 2 ) 


2-4 


2.7 


6.6±0.1±°; 2 


4-6.5 


5.0 


6.3±0.2+°; 2 


6.5-10 


7.8 


5.9 ± 0.2+°; 2 


10-15 


11.9 


5.5 ± 0.2t° i 


15-30 


19.7 


5.5 ± 0.3l°; 2 


30-80 


41.0 


4.9 ±0.618:1 



Table 3: The slope b resulting from a fit to the differential cross-section da/dt 
to an exponential form for the reaction 7*p — > p°p, for different Q 2 intervals. The 
first column gives the Q 2 bin, while the second column gives the Q 2 value at which 
the differential cross sections are quoted. The first uncertainty is statistical, the 
second systematic. 



Q 2 bin 


W bin 


Q 2 


W 




a(Yp - 


■+ p°p) 


(GeV 2 ) 


(GeV) 


(GeV 2 ) 


(GeV) 


(nb) 


stat. 


syst. 


2-3 


40- 


100 


2.4 


90 


647.1 


±8.7 


+28.4 
-41.7 


3-4 


40- 


100 


3.4 


90 


396.7 


±6.7 


+14.6 
-19.4 


4-5 


40- 


100 


4.4 


90 


247.8 


±5.8 


+8.9 
-12.6 


5-7 


40- 


120 


5.8 


90 


140.3 


±2.6 


+3.9 
-5.9 


7-10 


40- 


140 


8.2 


90 


71.9 


±1.4 


+1.7 
-2.9 


10-15 


40- 


140 


12 


90 


29.73 


±0.68 


+0.75 
-1.14 


15-20 


40- 


140 


17 


90 


12.77 


±0.50 


+0.27 
-0.42 


20-30 


40- 


140 


24 


90 


6.03 


±0.31 


+0.37 
-0.13 


30-50 


40- 


140 


37 


90 


1.88 


±0.16 


+0.07 
-0.15 


50-80 


40- 


140 


60 


90 


0.36 


±0.07 


+0.04 
-0.03 


80-160 


40- 


140 


100 


90 


0.05 


±0.03 


+0.02 
-0.01 



Table 4: Cross-section measurements at Q 2 and W = 90 GeV averaged over 
the Q 2 and W intervals given in the table. The normalisation uncertainty due to 
luminosity (±2%) and proton-dissociative background (±4%) ^ s n °t included. 
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Q 2 bin 
(GeV 2 ) 


W bin 

(GeV) 


Q 2 

(GeV 2 ) 


W 
(GeV) 


(nb) 


stat. 


■+ p°p) 
syst. 


2- 


-3 


32-40 


2.4 


36.0 


451.9 


±15.1 


+25.5 
-43.6 


2- 
2- 


-3 
-3 


40-60 
60-80 


2.4 
2.4 


50.0 
70.0 


554.1 
599.9 


±11.5 
±13.9 


+31.6 
-39.2 
+28.5 
-38.5 


2- 
2- 


-3 
-3 


80-100 
100-120 


2.4 
2.4 


90.0 
110.0 


622.5 
690.1 


±17.3 
±30.3 


+33.8 
-43.2 
+40.9 
-66.9 


3- 
3- 


-5 
-5 


32-40 
40-60 


3.7 
3.7 


36.0 
50.0 


240.8 
277.5 


±8.0 
±5.9 


+9.5 
-15.5 
+12.2 
-15.3 


3- 


-5 


60-80 


3.7 


70.0 


303.7 


±7.3 


+11.1 
-14.4 


3- 


-5 


80-100 


3.7 


90.0 


344.6 


±9.4 


+10.4 
-17.2 


3- 


-5 


100-120 


3.7 


110.0 


404.7 


±15.5 


+15.2 
-22.5 


5- 


-7 


32-40 


6.0 


36.0 


88.5 


±5.1 


+6.0 
-4.1 


5- 


-7 


40-60 


6.0 


50.0 


104.9 


±3.6 


+3.6 
-6.9 


5- 


-7 


60-80 


6.0 


70.0 


113.6 


±4.1 


+6.0 
-3.9 


5- 


-7 


80-100 


6.0 


90.0 


127.6 


±4.9 


+4.0 
-5.8 


5- 


-7 


100-120 


6.0 


110.0 


144.0 


±6.1 


+8.6 
-8.4 


7- 


-10 


40-60 


8.3 


50.0 


52.3 


±1.9 


+ 1.7 
-2.7 


7- 


-10 


60-80 


8.3 


70.0 


61.7 


±2.4 


+2.1 
-2.9 


7- 


-10 


80-100 


8.3 


90.0 


70.1 


±2.9 


+2.0 
-3.3 


7- 


-10 


100-120 


8.3 


110.0 


75.2 


±3.4 


+3.1 
-3.0 


7- 


-10 


120-140 


8.3 


130.0 


87.5 


±4.7 


+2.5 
-4.1 


10- 


-22 


40-60 


13.5 


50.0 


16.4 


±0.6 


+0.6 
-0.7 


10- 


-22 


60-80 


13.5 


70.0 


20.2 


±0.8 


+0.8 
-0.7 


10- 


-22 


80-100 


13.5 


90.0 


21.9 


±0.9 


+0.7 
-0.9 


10- 


-22 


100-120 


13.5 


110.0 


24.3 


±1.1 


+0.9 
-1.2 


10- 


-22 


120-140 


13.5 


130.0 


27.7 


±1.4 


+0.9 
-1.0 


10- 


-22 


140-160 


13.5 


150.0 


30.7 


±2.3 


+1.2 
-1.1 


22- 


-80 


40-60 


32.0 


50.0 


1.5 


±0.2 


+0.2 
-0.1 


22- 


-80 


60-80 


32.0 


70.0 


2.3 


±0.2 


+0.1 
-0.1 


22- 


-80 


80-100 


32.0 


90.0 


2.6 


±0.3 


+0.3 
-0.2 


22- 


-80 


100-120 


32.0 


110.0 


3.6 


±0.4 


+0.1 
-0.3 


22- 


-80 


120-140 


32.0 


130.0 


4.0 


±0.5 


+0.2 
-0.4 


22- 


-80 


140-160 


32.0 


150.0 


4.2 


±0.6 


+0.2 
-0.4 


22- 


-80 


160-180 


32.0 


170.0 


3.6 


±0.7 


+0.3 
-0.3 



Table 5: Cross-sections values obtained at Q 2 and W as a result of averaging 
over bins of the Q 2 and W intervals given in the table. The normalisation uncer- 
tainty due to luminosity (±2%) and proton-dissociative background (±4%)> are n °t 
included. 
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Q 2 bin 
(GeV 2 ) 


Q 2 

(GeV 2 ) 


5 stat. syst. 


2- 3 

3- 5 
5-7 
7-10 

10-22 
22-80 


2.4 
3.7 
6.0 
8.3 
13.5 
32.0 


0.321 ±0.035 toitl 
0.412 ±0.036 j£o35 
0.400 ±0.052 
0.503 ±0.057 
0.529 ±0.051 ±g:ggg 
0.834 ±0.118 !Hi2 



Table 6: TTie mfote o/ 5 obtained from fitting H^ 5 . TTie /irsi column 

gives the Q 2 bin, while the second column gives the Q 2 value at which the cross 
section was quoted. 



Q 2 bin (GeV 2 ) 


Q 2 (GeV 2 ) 


W bin (GeV) 


r 04 
' 00 


R = a L /a T 


2-3 


2.4 


32 - 120 


0.60 ± O.Oll&M 


1 ccn+0-05 +0.20 
1,JU -0.05 -0.15 


3-5 


3.7 


32 - 120 


0.68 ±0.01+^ 


9 1 n +0.08 +0.18 
z - iu -0.08 -0.14 


5-7 


5.9 


40 - 140 


0.73 ±0.01+^ 


o 7n +0.14 +0.26 
z - ' u -0.13 -0.28 


7-10 


8.3 


40 - 140 


0.76 ± 0.01+{j;{£ 


o 9 n +0.20 +0.25 
°- zu -0.18 -0.27 


10-15 


12.0 


40 - 140 


0.78 ±0.0llg;gl 


o cr n + - 26 + - 30 

O.OU_ 2 4 _o.26 


15-30 


19.5 


40 - 140 


0.82 ± 0.02+°;^ 


a fin+ - 54 + - 48 
^.ou_ 45 _ 44 


30-100 


40.5 


40 - 160 


0.86 ± 0.04+°;^ 


c m+ 2 - 75 + 2 - 15 

u - iu -1.56 -0.85 



Table 7: The spin matrix element and the ratio of cross sections for longitudi- 
nally and transversely polarised photons, R = olI^t, as a function of Q 2 , averaged 
over the Q 2 and W bins given in the table. The first uncertainty is statistical, the 
second systematic. 
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Q 2 bin (GeV 2 ) 


Q 2 (GeV 2 ) 


W bin (GeV) 


W (GeV) 


r 04 
' 00 


R = a L /a T 


2-3 


2.4 


32-55 


43 


0.60 ± O.Ollara 


1 crn+0.06 +0.21 
i - JU -0.06 -0.15 


2-3 


2.4 


55-75 


65 


0.60 ± 0.01±°;°3 


i cn+ - 06 + - 35 

1 - du -0.06 -0.17 


2-3 


2.4 


75-110 


91 


0.59 ± 0.0llg;{£ 


1 /iq+0.06 +0.23 
i -^°-0.06 -0.23 


3-7 


4.2 


32-60 


45 


0.70 ± O.Oll^ 


9 OO+0.09 +0.13 
z -°°-0.09 -0.09 


3-7 


4.2 


60-80 


70 


0.69±0.0ll°m 


o 90+O.I2 +0.24 
z - zo -0.11 -0.10 


3-7 


4.2 


80-120 


99 


0.69 ±0.01^} 


9 90+O.IO +0.14 
z - zo -0.09 -0.09 


7-12 


8.8 


40-70 


55 


0.74 ±0.0llg;^ 


9 oc+0.25 +0.23 
Z.OO„ 22 -0.26 


7-12 


8.8 


70-100 


85 


0.76 ± 0.02lg;{£ 


1 7+0.38 +0.19 

01 ' -0.32 -0.28 


7-12 


8.8 


100-140 


120 


0.76±0.02l°;°2 


1 7+0.38 +0.23 

01 ' -0.32 -0.26 


12-50 


18.0 


40-70 


55 


0.84 ± 0.03tg$ 


r nc+1.16 +0.54 
J - ZJ -0.84 -0.34 


12-50 


18.0 


70-100 


85 


0.82 ± 0.03lg;{£ 


A cc+0.94 +0.47 
^• dd -0.70 -0.43 


12-50 


18.0 


100-160 


130 


0.83 ± 0.02+°;° 2 


A qq+0.87 +0.64 
^• OS -0.67 -0.39 



Table 8: The spin matrix element and the ratio of cross sections for longi- 
tudinally and transversely polarised photons, R = o L ja T , as a function of W for 
different values of Q 2 , averaged over the Q 2 and W bins given in the table. The 
first uncertainty is statistical, the second systematic. 
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O 2 bin (GeV 2 ) 


O 2 (GeV 2 ) 


W bin (GeV) 

V V Ulll I V.I Vj V J 




t\ (GeV 2 ) 


r 04 
' 00 


r? = CT j / fjrn 


2-5 


3.0 


32 - 120 


0.04 


0.62 ±0.01^ 


i ^q+0.07 +0.15 
1,uo -0.06 -0.13 


2-5 


3.0 


32 - 120 


0.14 


0.62 iO.Oll&ra 


i co+0.09 +0.10 
±.UO_qq 9 _Q.i9 


2-5 


3.0 


32 - 120 


0.27 


0.63 ±0.0lljj;{£ 


1 7H+0-H +0.24 
L - 1 u -0.11 -0.14 


2-5 


3.0 


32 - 120 


0.45 


0.64 ± 0.02+°;° 2 


i 70+O.I4 +0.16 

1 . 1 0_n i o „n 91 

— u.lj — Ij . £1 ± 


2-5 


3.0 


32 - 120 


0.76 


0.63 ± 0.03ljj;{£ 


1 7n +0.26 +0.63 
L - ' u -0.22 -0.32 


5-50 


10.0 


40 - 160 


0.04 


0.74 ±0.0llg;gl 


9 O/1+0.18 +0.16 
z -°^-0.17 -0.15 


5-50 


10.0 


40 - 160 


0.15 


0.75 ±0.01+^2 


nn+ - 26 + - 17 

O.UU_ 23 _ 30 


5-50 


10.0 


40 - 160 


0.27 


0.74 ± 0.02+°;° 2 


9 O/1+0.26 +0.32 
z -°^-0.24 -0.51 


5-50 


10.0 


40 - 160 


0.45 


0.72 ± 0.02+°;°^ 


9 p-7+0.29 +0.41 
' -0.25 -0.22 


5-50 


10.0 


40 - 160 


0.76 


0.73 ± 0.04+°;^ 


9 7f|+ - 56 + - 45 
z - ' u -0.43 -0.57 



Table 9: The spin matrix element and the ratio of cross sections for longi- 
tudinally and transversely polarised photons, R = olIo~t, as a function of \t\ for 
two values of Q 2 , averaged over the Q 2 and W bins given in the table. The first 
uncertainty is statistical, the second systematic. 



Q 2 bin (GeV 2 ) 


Q 2 (GeV 2 ) 


|*| (GeV 2 ) 


a P (t) 


2-5 


3 


0.04 


1.104 ±0.0111°;°}° 


2-5 


3 


0.14 


1.099 ±0.0141°;°^ 


2-5 


3 


0.28 


1.048 ±0.016iH?4 


2-5 


3 


0.57 


1.013 ±0.0211^17 


5-50 


10 


0.04 


1.149 ±0.0121°;°^ 


5-50 


10 


0.16 


1.134 ±0.0141°,;°$ 


5-50 


10 


0.35 


1.104 ±0.017l°;°} 2 


5-50 


10 


0.68 


1.085 ± 0.028l°;°^ 2 



Table 10: The values of the effective Pomeron trajectory a P (t) as a function of 
\t\, for two Q 2 values. The first uncertainty is statistical, the second systematic. 
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Q 2 bin (GeV 2 ) 


Q 2 (GeV 2 ) 


a F {0) 


c4(GeV" 2 ) 


2-5 
5-50 


3 

10 


1.113 ±0.010+°;°^ 
1.152 ±0.0111^6 


0.185 ± 0.042+°;° 22 
0.114 ±0.0431Hm 



Table 11: The values of the effective Pomeron trajectory intercept a P (0) and slope 
a'p, for two Q 2 values. The first uncertainty is statistical, the second systematic. 



Q 2 (GeV 2 ) 


W (GeV) 


b (GeV" 2 ) 


3.5 


38 


6.3 ± 0.2±° i 


3.5 


57 


6.3±0.ll°j 


3.5 


82 


6.6±0.2+°; 2 


3.5 


107 


6.9 ± 0.2lgj 


3.5 


134 


7.0 ± 0.3ig;| 


11 


38 


5.8 ± O.Stoi 


11 


57 


5.8 ± 0.2+°; 2 


11 


82 


5.7±0.2t°; 2 


11 


107 


5.9 ± 0.2+^ 


11 


134 


6.1 ±0.2+°;^ 



Table 12: The slope b resulting from a fit of the differential cross section da/dt 
for the reaction 7*p — > p°p to an exponential form, for different W values, for two 
Q 2 values. The first uncertainty is statistical, the second systematic. 
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Figure 1: Comparison between the data and the Zeusvm MC distributions for (a) 
W, (b) Q 2 , (c) \t\, (d) cos6 h , (e) <p h and (f) & h for events with 0.65 < M n7r < 1.1 
GeV and \t\ < 1.0 GeV 2 . The MC distributions are normalised to the data. 
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Figure 2: Comparison between the data and the Zeusvm MC distributions for 
the transverse momentum, p?, of ir + and ir~ particles, for different ranges of Q 2 , 
as indicated in the figure. The events are selected to be within 0.65 < M n7T < 1.1 
GeV and \t\ < 1.0 GeV 2 . The MC distributions are normalised to the data. 
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Figure 3: (a) The energy distribution in the FPC. The data (full dots) are 
compared to the expectations from the Pythia MC, normalised to the data, (b) The 
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from the Pythia MC, normalised to the data for xl < 0.95. The extracted fraction 
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(open circles), as a function of (c) Q 2 , (d) W and (e) \t\. All events were selected 
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Figure 5: The 7r + 7r~ acceptance- corrected invariant-mass distribution, for different 
Q 2 intervals, with mean values as indicated in the figure. The lines are defined in 
the caption of Fig. [7J 
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Figure 6: The 7r + 7r" acceptance- corrected invariant-mass distribution, for different 
t intervals, with mean values as indicated in the figure. The lines are defined in the 
caption of Fig. [7J 
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Figure 7: The 15 density-matrix elements obtained from a fit to the data (dots), 
as a function of Q 2 . Also shown in the figure are results from an earlier measure- 
ment [50] (open circles). The inner error bars indicate the statistical uncertainty, 
the outer error bars represent the statistical and systematic uncertainty added in 
quadrature. The dotted line at zero is the expectation from SCHC when relevant. 
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Figure 8: The acceptance- corrected cos 9 'h distribution, for different Q 2 intervals, 
with mean values indicated in the figure. The line represent the fit to the data of 
Eg. 
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Figure 9: The differential cross-section da / d\t\ as a function of\t\ for^p — > p°p, 
for fixed values of Q 2 , as indicated in the figure. The line represents an exponential 
fit to the data. The inner error bars indicate the statistical uncertainty, the outer 
error bars represent the statistical and systematic uncertainty added in quadrature. 
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Figure 10: The value of the slope b from a fit of the form do~/d\t\ oc e~ fe '*' for 
exclusive p° electroproduction, as a function of Q 2 . Also shown are values of b 
obtained previously at lower Q 2 values [10, 53]. The inner error bars indicate the 
statistical uncertainty, the outer error bars represent the statistical and systematic 
uncertainty added in quadrature. 
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Figure 11: A compilation of the value of the slope b from a fit of the form 
da/d\t\ oc e~ 6 '*' for exclusive vector-meson electroproduction, as a function of Q 2 + 
M 2 . Also included is the DVCS result. The inner error bars indicate the statistical 
uncertainty, the outer error bars represent the statistical and systematic uncertainty 
added in quadrature. 
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Figure 12: The Q 2 dependence of the cross section for exclusive p° electropro- 
duction, at a 7*p centre- of -mass energy W =90 GeV. The ZEUS 1994 [53] and 
the ZEUS 1995 [10] data points have been extrapolated to W = 90 GeV using the 
parameterisations reported in the respective publications. The inner error bars in- 
dicate the statistical uncertainty, the outer error bars represent the statistical and 
systematic uncertainty added in quadrature. 
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Figure 13: The W dependence of the cross section for exclusive p° electropro- 
duction, for different Q 2 values, as indicated in the figure. The inner error bars 
indicate the statistical uncertainty, the outer error bars represent the statistical and 
systematic uncertainty added in quadrature. The lines are the result of a fit of the 
form W s to the data. 
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Figure 14: The value of 5 from a fit of the form W 5 for exclusive p° electropro- 
duction , as a function of Q 2 . Also shown are values of 5 obtained previously at 
lower Q 2 values [10, 53]. The inner error bars indicate the statistical uncertainty, 
the outer error bars represent the statistical and systematic uncertainty added in 
quadrature. 
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Figure 15: Comparison of the HI (squares) and ZEUS (dots) measurements of the 
W dependence of a' 1 * v ^' p p , for different Q 2 values, as indicated in the figure. The 
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Figure 16: A compilation of the value of 5 from a fit of the form W s for exclusive 
vector-meson electroproduction, as a function of Q 2 + M 2 . It includes also the 
DVCS results. The inner error bars indicate the statistical uncertainty, the outer 
error bars represent the statistical and systematic uncertainty added in quadrature. 
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Figure 17: The ratio as a function of Q 2 for W = 90 GeV. Also included are 
values o/ Tqq from previous measurements at lower Q 2 values [10, 53]. The inner 
error bars indicate the statistical uncertainty, the outer error bars represent the 
statistical and systematic uncertainty added in quadrature. 
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Figure 18: Comparison of the HI (squares) and ZEUS (dots)measurements of 
R as a function of Q 2 . The HI data are at W = 75 GeV and those of ZEUS at 
W = 90 GeV. Also included are measurements performed previously at lower Q 2 
values [10, 53]. The inner error bars indicate the statistical uncertainty, the outer 
error bars represent the statistical and systematic uncertainty added in quadrature. 
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Figure 19: The ratio R as a function of M n7T! for W = 80 GeV, and for two 
values of Q 2 , as indicated in the figure. The inner error bars indicate the statistical 
uncertainty, the outer error bars represent the statistical and systematic uncertainty 
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Figure 20: The ratio as a function ofW for different values ofQ 2 , as indicated 
in the figure. The inner error bars indicate the statistical uncertainty, the outer 
error bars represent the statistical and systematic uncertainty added in quadrature. 
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Figure 21: TTie ratto Tqq as a function of\t\ for different values ofQ 2 , as indicated 
in the figure. The inner error bars indicate the statistical uncertainty, the outer 
error bars represent the statistical and systematic uncertainty added in quadrature. 
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Figure 22: The effective Pomeron trajectory a P (t) as a function oft, for two 
values of Q 2 , with average values indicated in the figure. The inner error bars 
indicate the statistical uncertainty, the outer error bars represent the statistical and 
systematic uncertainty added in quadrature. 
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Figure 23: The parameters of the effective Pomeron trajectory in exclusive p° 
electroproduction, (a) a P (0) and (b) a' p , as a function of Q 2 . The inner error bars 
indicate the statistical uncertainty, the outer error bars represent the statistical and 
systematic uncertainty added in quadrature. The band in (a) and the dashed line 
in (b) are at the values of the parameters of the soft Pomeron [19, 20]. 
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Figure 24: The b slope as a function ofW for two ranges ofQ 2 , with average values 
as indicated in the figure. The inner error bars indicate the statistical uncertainty, 
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Figure 25: The Q 2 dependence of the 7*p — > p°p cross section at W =90 GeV. 
The same data are plotted in (a) and (b), compared to different models, as described 
in the text. The predictions are plotted in the range as provided by the authors. 
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Figure 26: The W dependence of the 7*p — > p°p cross section for different 
values of Q 2 , as indicated in the figure. The same data are plotted in (a) and (b), 
compared to different models, as described in the text. The predictions are plotted 
in the range as provided by the authors. 
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Figure 27: The value of 5 from a fit of the form a ~ W 5 for the reaction 
7*p — > p°p, as a function of Q 2 . The lines are the predictions of models as denoted 
in the figure (see text). 
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Figure 28: 

reaction ^*p — > p°p, as a function of Q 
as denoted in the figure (see text). 
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Figure 29: The ratio as a function of Q 2 compared to the predictions of 
models as denoted in the figure (see text). 
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Figure 30: The ratio Tqq as a function ofW for different values of Q 2 compared 
to the predictions of models as indicated in the figure (see text). 
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Figure 31: The ratio Tqq as a function ofW for different values of Q 2 compared 
to the predictions of models as indicated in the figure (see text). 
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